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Abstract. We show with a short Chandra HRC-S exposure that the powerful X-ray source coincident with the nucleus of 
M33 is unresolved at the highest spatial resolution available to date (0'.'4). The source flux is variable, doubling during the 
5 ks exposure. The combination of properties exhibited by M33 X-8 establishes the source as the nearest example of the ultra- 
luminous X-ray sources that have been uncovered in other nearby galaxies. On short timescales, we set limits of 9% rm.s. 
variability for pulsations in the 0.01-1000 Hz range. 
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1. Introduction 



The nearby galaxy M33 (d ^ 795 kpc; Ivan den Berghl 
Il99ll) hosts the most luminous steady X-ray source in 
the Local Group: M33 X-8 has a 1-10 keV lumi- 
nosity of 1.2-10^^ erg s"'. The source has been de- 
tected at this l evel ever since the first X -ray observa- 
tions of M33 (Long. D'Odorico. Charles. & Dopita. 1981; 



(iLong^ D'Odorico. Charles. & Dopita, 
sl ll983l) . 



M33 

iMarkert & RalU 

Chandra observations established that M33 X-8 is coinci- 
dent with the optical position of th e nucleus to the (y.'6 uncer- 
tainty of the astrometric solution (iDubus & Rutledgel 12002). 
However, the 1500 M© upper limit on the mass of a cen- 
tral black hole placed by the measured velocity dispersion 
(JGebhardt et all 1200 ll) an d the lack of activity at other wave- 
lengths JLong et alJ200 2 and references therein) rule out a low 
luminosity active galactic nucleus. 

Ultra-luminous X-ray sources (ULXs) are defined as point 
sources having X-ray luminosities in excess of the Eddington 
limit for a 10 M© object (10^'^ erg s"'; Fabbiano & Whitj 
I2OO3I) . ULXs could be due to beamed emission from stel- 
lar mass-sized neutron stars (NS) / black holes (BH) or truly 
isotropic emission from 100-1000 Mq intermediate mass BHs 
( JKingetallbOOll) . If M33 X-8 is a single object, then it is the 
nearest example of a ULX. The detection of a 106 day ~20% 
modulation in t he X-ray flux from M33 X-8 in a 5 year span 
of ROSAT data. JDubuset ailll997l) and flie resemblance of the 
spectrum to that of black hole binaries ( iLa Parola et all 120031) 
support this picture. 

ROS AT images show M33 X-8 to be point-like at 5" res- 
olution (Schulman & Bregmanl Il995l) . But diffuse emission 
from the compact optical nucleus (3" wide, see Fig. 3 in 



iKormendv & McClurel 19931) or contributions from close, faint 
sources would not have been resolved using any pre-Chandra 
instrument. With its superior angular resolution, Chandra has 
been able to reduce the possibility that ULXs in nearby galax- 
ies are chance superpositions of lower luminosity sources. 
Unfortunately, in the case of M33 X-8, the existing Chandra 
observations of M33 X-8 had the nuclear source either far off- 
axis or heavily piled-up in ACIS CCD. As a result, the Chandra 
data could not be used to study the radial extent below a few 
arc seconds. 

Here we report a new observation of M33 X-8 with 
Chandra, an on axis observation obtained with the High 
Resolution Camera (HRC) on board Chandra that fully ex- 
ploits the spatial resolution of Chandra in an attempt to resolve 
the nuclear source in M33. The observation, spatial and timing 
analysis are described in §2-4. The results are discussed in §5. 

2. Observation 

M33 was observed with the ChandrafHRC-S instrument 
{Murray etal. , 1998 ) on July 29, 2003 from 19:53:12 UT for 
a total exposure time of 5355 s. The HRC is unaffected by 
pile-up and provides the best achievable spatial resolution (0'.'4 
FWHM) to date. We decided to use the HRC-S which has 
the capability to study the poorly known timing properties of 
M33 X-8 at frequencies above 1 Hz. We performed the analy- 
sis with ClAO version 3.0.1'. 

There is only one point source, M33 X-8, with signal-to- 
noise greater than 5 in the field-of-view. In the 6x8' cen- 
tral region surrounding M33 X-8 where the background is 



.http : //asc . harvard . edu/ciao 



G. Dubus et al.: High resolution Chandra X-ray imaging of the nucleus of M33 




2000 3000 

Time (s) 



100 



o 
sx 

>, 
Si 

CD 



10 




1 
10 



10- 



I I I L 

10-2 10-1 10" 101 
frequency (Hz) 



102 



103 



Fig.l. HRC-S lightcurve of M33 X-8 in 100s bins, showing 
significant variability on timescales of ~ 1 ks. 

spatially uniform, we estimated the average background rate 
per 2" circular region, and found that 28 counts would be 
needed for a 3cr detection (taking into account the 13750 sep- 
arate trials). This translates into a 0.6-9 keV luminosity up- 
per limit for point source detection in this observation of < 
9 10"'^(t//795 kpc)- erg s"' (assuming a photon power-law of 
a-2, Nh = 1.9 X 10^'cm"^). There are 43 so urces within 8' of 
M33X-8 in the X-ray ROSAT catalogue of 'Haberl & Pietsch' 



i2001l) but the brightest of these (excluding X-5 which is at the 
very edge of the field-of-view) have count rates a factor 5 1 % 
that of M33 X-8. Hence, it is not surprising that we failed to 
detect any other sources in this short exposure. 

3. Timing Analysis 

We extracted 3756 counts from a 16 pixel radius region around 
M33 X-8, giving an average of 0.70+0.01 counts s"'. Using 
an annulus with inner and outer radii of 20 and 30 pixels, we 
estimate the background contribution to be only 6 counts (and 
hence negligible). Surprisingly, the count rate from M33 X- 
8 varied significantly during our observation, increasing from 
0.36+0.03 to 0.92±0.04 counts s ', on a timescale of 2400 s 
(Fig. Q. Previous results have consistently shown that the X- 
ray spectrum of M33 X-8 is well fitted by a ^r=3.5 keV ther- 
mal bremsstrahlung spectrum wit h A^;:^ == 1.9 x 10^' cm~^ 
("Takano et al.', '1994.'; 'Parm ar et all |200l'; 'Pubus & Rut ledgeL 
12002: LaParolaet al.. 20031) . Assuming this spectrum (for 
which 1 HRC-S count corresponds to 2.6- 10 " erg cm"^ of 
unabsorbed flux in the 0.6-9 keV band), the maximum count 
rate is equivalent to Lx - 1.8-10"'' id 1195 kpc)^ erg s"'. 

We applied barycenter corrections to the photon arrival 
times and computed their power density spectrum (PDS), 
with a Nyquist frequency of 1000 Hz. We find no excess 



Fig. 2. Power Density Spectrum of M33 X-8 during the 
53 55 s Chandra HRC -S observation, normalized according 
to iLeahv et al.l ( 119831) . Data have been rebinned logarith- 
mically. The solid line shows the best-fit power-law (slope 
a=2. 65+0.15) + constant noise (taken equal to the average at 
high frequencies, 1.98). 

power above 0.01 Hz, with a 99.9% upper limit of <8.8% 
root-mean-squared (rm.s.) variability from a sinusoidal signal. 
Logarithmically rebinning the PDS to 100 bins, we do not find 
any excess from broad timing features above >0.01 Hz. The 
90% upper limit to a quasi-periodic oscillation (QPO) is 24% 
r.m.s. where we have modeled the QPO as a 0.25 Hz FWHM 
gaussian at 6 Hz. There is significant power at low frequen- 
cies with a measured 30+6% r.m.s. variability below 0.01 Hz. 
Further rebinning to 20 frequency bins (Fig.|2j, we find that the 
data are well described (^^=0.7, 18 degrees of freedom) by a fit 
to a constant plus a power-law low frequency noise (LFN) with 
a steep slope q'=2.65+0.15 (90% confidence error). The con- 
stant is taken equal to the high frequency average, 1 .98 (close 
to the expected level from a pure Poisson flux). However, the 
LFN power must turn over between 0.02 and 0.19 mHz in or- 
der for it not to diverge (i.e. exceed 100%). To constrain any 
additional high frequency, flatter timing feature, we fixed the 
power-law slope at -2.65 and added a second component with 
a - \. The 90% confidence upper limit to the rm.s. variabil- 
ity below 10 Hz is <3.2%, corresponding to a change in the 
minimum ;t'^ of 4.61 (iLampton. Margon. & Bowver 111 9761) . 

In spite of the lack of photon energy sensitivity in the HRC 
detectors, we can still search for spectral variability by exploit- 
ing the fact that the photon energy to pulse-invariant (PI) en- 
ergy channel calibration is stable during an observation. We 
therefore investigated the data in two ways: firstly by compar- 
ing counts in the first 2000 s interval (cf. Fig.Q) with those de- 
tected afterwards (a total of 907 and 2829 counts respectively); 
secondly by dividing in time so as to produce an equal number 
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of counts in each group (which occurs 3172 s into the obser- 
vation). The distributions of the PI channels between groups 
are statistically identical (respectively a probability of 0.30 and 
0.19 of being drawn from the same distribution for a two-sided 
Kolmogorov-Smirnov test). Hence, within the constraints of 
the HRC, there is no evidence for spectral variability during^ 
our observation. ^ 

o 

4. Spatial Analysis 5 

Our main purpose in undertaking this observation was to deter- g 
mine whether the nuclear source in M33 is point-like. We there- 
fore used ChaRT and Marx to simulate how the PSF of M33 X- 
8 would appear on the HRC-S detector if it were point-like. In 
carrying out this simulation (of 28861 counts in total), we as- 
sumed the same spectrum as in §3. This simulation accounts 
for all known telescope and detector effects. It also shows why 
the spectrum is important, as the ChaRT PSF is significantly 
broader than a monochromatic PSF (at 1 keV) calculated with 
mkpsf. 

We logarithmically re-binned the data into a 1-D radial dis- 
tribution (between r = and 40") for comparison with a ChaRT 
simulation (Fig.|3j- We determined the source centre by rebin- 
ning the PSF and performing a radial fit to the data, minimis- 
ing x^ against the 2-D offset. The best-fit to the PSF model 



has ^y= 1.72 (18 degrees of freedom), and the simulated PSF 
is a good description of the data. We find no evidence for any 
additional component to the emission profile. 

Extended emission associated with the nucleus is not sta- 
tistically required by the data. To place an upper limit on its 
contribution, we convolved the PSF with the radi al profile de- 
rived from Hubble Space Telescope visible data JDubus et all 
11999 *). The 90% confidence upper limit is <3.7% of the average 
flux from M33 X-8 during the observation or 4.7 ■ 10^^ erg s"' . 

We also examined the residuals obtained by subtracting 
the best-fit PSF from the data using an image bin size of 3 
pixels (0'.'4), equal to the resolution of the HRC-S. There is 
no apparent structure in the residuals, which have amplitudes 
smaller than <7% of the peak flux. We derive an upper limit of 
9. 1 -10^^ erg s"' to the contribution from a point source between 
O'.'4-l" from the nucleus. 

5. Discussion 

[Dubus & Rutledge ( 2002) established to a high degree of pre- 
cision that M33 X-8 is coincident with the nucleus. The present 
Chandra HRC-S observation shows that M33 X-8 is also un- 
resolved at 0'.'4 resolution, improvin g on earli er ROSAT HRI 
upper limits by a factor ~ 10 ( Schulman & Bregman. 1995) . 
This rules out contributions from nearby sources but does not 
rule out mul tiple sources within the compact nucleus core 
toubus et all 1 1999) . However, the substantial (50%) variabil- 
ity we observe during this short 5 ks observation shows that a 
single source dominates. We conclude M33 X-8 is most likely 
a single source and, as has been claimed previously, the nearest 
example of a ULX source. 

A possible key to uncovering the nature of this source is 
its variability. The short timescale doubling that we observe 



1000 



100 r 



10 r 



1 r 

4 r" 

2 r 
X Or' 

-2 r 
-4 _ 
0.1 






1.0 10.0 

arcsec 



Fig. 3. Top Panel: The observed radial profile of M33 X-8, 
showing the number of counts (with Gaussian errors) per unit 
area in annular bins. The dashed line is the best fit model 
(from ChaRT), and includes a constant background component. 
M33 X-8 dominates at radii <4", while the background dom- 
inates at larger radii. Bottom Panel: The ;^f=(data-model)/c^ 
residuals. We attribute the scatter between radii of 5 and 40"as 
due to unresolved point sources. 



is unprecedented: within single pointed observations M33 X- 
8 is usually seen to be constant. Variability at a level of 
10% has been obs erved within h ours to days in only two 
ROSAT visits (Dub us et all Il997l) . More recently, 10% vari- 
ability on a timescale of 5000 s was seen in a continuous 93 ks 
Chandra off-axis observation (iLa Parola et al...2003: see also 
IPeres. Reale. Collura. & Fabbianolll98^ for an earlier report). 
Although the timescale of the variation we observe is similar, 
the amplitude in our observation is rnuch la rger than the 10% 
variation observed bv lLa Parola et all ( l2003l) . 

Using the doubling timescale of 2400 s, we can set an upper 
limit to the emission region of 7-10'"' cm which hints at a com- 
pact object. Recently, it has been proposed that the break fre- 
quency at which the PDS of active galactic nuclei and X-ray bi- 
naries switch fro m a -2 to a -1 power-la w slope depends on the 
black hole mass JMcHardv et all 120041) . Using the fact that the 
PDS of M33 X-8 must turn over between 0.02 and 0.19 mHz, 
and interpr eting this break as the high frequency turnover (e.g. 
NGC 405 1 . iMcHardv et alJ2004h . we find the black hole mass 
would be in the range 10^"^ M^. This is inconsistent with ve- 
locity dispersion upper limits (Gebhardt et al.. 2001) . Hence, it 
is unlikely the variability is asso ciated with the high fre quency 
noise phenomenon described by McHa rdv et all 12004). Much 
longer exposures are needed to establish the timing properties 
of M33 X-8 and investigate its nature. 

The X-ray spectrum, long timescale variability and lack of 
an optical counterpart suggest an X-ray binary. The ~1 hour 
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doubling, but lack of fast variability, are consistent with a low 
mass X-ray binary in the hi g h or v ery high state. On long 
timescales, iMarkert & RallisI (Il983h observed a 50% count 
rate drop between two Einstein HRI observations separated 
by six months. Later ROSAT ohservations confirmed variation s 
of k, 20% on a timescale of 106 days ("Pub us et all Il997l) . 
iLa Parola et alJ ( I2003.) have emphasized the similarity between 
M33 X-8 and LMC X-3. LMC X-3 has long term varia tions 
on a 99 day (or 199 day) timescale (Wil ms et all l200lh . We 
note that the present observation would correspond to phase 
cp = 0.87 in Fig. 3 of Dubus et al. ( 1997), very close to the 
minimum of the cycle where the folded ^0&4rdata shows the 
most scatter. 

Chandra results indicate that ULXs are more fre- 
quently found in assoc iation with recent star formation 
JFabbiano & WhitelEoOSi) . There is evidenc e for a 40 Myr old 
starbu rst in the nucleus of M33 (lb'Connelll ll983HLong et all 
I2OO2I) . The presence of M33 X-8 could be linked to this episode 
of star formation, suggesting perhaps a source comprised of 
a high mass companion orbiting a stellar-mass black hole. At 
present, the source in the nucleus of M33 appears to have all 
of the characteristics of a canonical ULX. As such, a better un- 
derstanding of the nearest ULX should remain a priority. 
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